Extension of Frohlich's method to 4-fermion interactions 
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Higher order terms of the transformed electron-phonon Hamiltonian H e - P h, obtained by perform- 
ing the Frohlich's transformation, are investigated. The influence of terms discarded by Frohlich 
(in particular those proportional to the third power of electron-phonon coupling) on the effective 
Hamiltonian is examined. To this end a second Frohlich-type transformation is performed, which 
yields, among others, an effective 4-electron interaction. This interaction is reduced to a form ad- 
mitting solution of thermodynamics. The form of the coupling of the 4-electron interaction is found. 
By applying standard approximations, it is shown that this interaction is attractive with interaction 
coupling given by — £)£ /w£ , where Dk is electron-phonon coupling, lo^ is phonon energy and kF 
is Fermi momentum. 

The form of higher order terms of the original Frohlich-transformed H e - P h are also found, up to 
terms proportional to the 6-th power of the coupling, that is up to those, which yield the effective 
4-electron interactions. 
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I. INTRODUCTION 

The BCS theorji proved extremely effective as a theory of superconductivity. The key idea of this theory — 
the attractive interaction which binds electrons into Cooper pairs — has its essential origins in earlier conclusions drawn 
Frohlich^, who performed a unitary transformation U of the electron-phonon Hamiltonian iJ e _ p /j. U was adjusted so 
as to eliminate the electron-phonon interaction as far as possible and replace it by an effective interaction V e f / between 
electrons dressed in the phonon field. V e // proved to be attractive for 1-electron energies close to Ef- A reduced form 
of V e f f was subsequently used by Bardeen, Cooper and Schrieffer in their theoryi. 

It is worth emphasizing that Frohlich's transformation is not strictly unitary, because the less significant terms of 
the resulting expansion of I) ' H e - v hXJ~ x were discarded. Correctness of the remaining terms, included into the Hamil- 
tonian of a superconductor, was confirmed by the success of the BCS theory. 

Unfortunately, BCS theory proved incapable of explaining superconductivity in type-II superconductors, heavy 
fermions and high-T c superconductors (HTSC). The search for an alternative theory of superconductivity proceeds in 
various directions and one of them exploits the idea of extending BCS theory by adding to the BCS Hamiltonian Hbcs 
further interactions. Rickayzen^ suggested incorporation of a 4-fermion interaction, motivating such choice by analo- 
gies between theory of superconductivity and nuclear theory, where such interactions had been considered. This 
idea was also mentioned by Volovik**, but remained in the realm of theoretical concepts until the late nineties, when 
Mackowiak and one of authors (PT)Wi 8 iV°i 1:L (MT) proposed a Hamiltonian H M t = H BC s + W + V M T, where 

Hbcs — T + V B cs = ^^ knkCT - I^T* ^kk'ak+ a -k- a -k'-flk'+, (1) 

ko- kk' 

W = J^7k^k+"k-, (2) 

k 



Vmt = |A| 2^ gkk'a k - a k+ a -k- a -k+ a -k'+a-k'-ak'+ak'-, (3) 

k,k' 

|A| in Eqs. denotes the system's volume, whereas guk>, Gkk' are bounded functions and 7k, ffkk', Gkk' are 

invariant under time reversal k — > — k or k' — > — k'. 

Considerations which led to the introduction of W were founded on the analysis of the HTSC normal state. 
There are grounds to believe that this is not a normal Fermi liquid state. The interaction W, first added to Hbcs 
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by Czerwonkc 12,13 , guaranteed normal-state behaviour characteristic of the so-called statistical spin liquid, considered 
earlier by Spalek and Wojci k 14 ' 15 . 

W is a 2-electron interaction. This raises the question, whether it can be obtained by a reduction procedure 
(different than the BCS one) of the interaction derived by Frohlich. It can be easily verified that this is impossible. 
More precisely, for the unique possible reduction of momenta, the coupling vanishes, meaning that the nature of W 
is not phononic. 

Introduction of the 4-fermion interaction Vmt was justified in Ref. |^ by its possible role as an attraction between 
pairs in HTSC, mediated by phonons or other quanta™. A conjecture put forward in Ref. |5| suggested also that Vmt 
could be expected to arise as one of the higher-order terms of Frohlich's expansion of UH e - p hU~ 1 . An alternative 
justification was given in Ref. |H where Vmt was viewed as a BCS-type interaction between quasi-particles of a free 
gas represented by W written in the form X)kTkC k Ck, with Ck = Ok+ak-- Both of these ideas essentially exploit 
the concept of phonon-type mediation of interactions. The significance of this mediation in HTSC was stressed 
by Wysokihskii&. 

The question of the form of higher-order terms of Frohlich's expansion is interesting itself, not only as providing 
a possible explanation of the MT extension, but, first of all, because these terms could throw some light on further 
possible extensions of Hbcs related to effects of electron-lattice interaction. Since Fermi-liquid theory will remain 
the foundation of our formalism, we shall focus our interest on effective electron interactions. 

It is worth noting that the possible presence of fcrmion quadruples in superconductors and superfiuids was considered 
in a number of papers. Schneider and Keller— measured the various characteristics of some cuprates and Chevrel- 
phase superconductors, especially concentrating on the relation between the critical temperature and zero temperature 
condensate density. They noticed that the experimental data for e.g., YBa 2 Cu3Og.602 point to similarities with 
the behaviour of a dilute Bose gas. As a result they suggested Bose condensation of weakly interacting fermion pairs 
as a mechanism of transition from normal to superconducting state. Bunkov et al— pointed to presence of fermion 
quadruples in 3 He. Their work was devoted to the problem of influence of spatial disorder on the order parameter in 
superfluid 3 He. By resorting to the work of Volovik^, they suggested that unusual spectra of 3 He in aerogel could be 
explained by a process in which impurities tend to destroy the anisotropic correlations of the order parameter, while 
correlations of higher symmetry can survive (e.g., four-particle correlations). Recently Schneider et al— discovered 
of half-/i/2e magnetic flux in SQUIDs fabricated of bicrystalline YBa2Cu307_a films. This situation corresponds to 
the presence of fermion quadruples in the system. Based on this observation, Aligia et al— investigated a model of 
interface between two superconductors, based on a one dimensional boson lattice model and proposed formation of 
quartets of electrons. 

In SectionlTTlhigher order terms of the expansion of Frohlich's transformation /7ff e _ p ^J7 _1 are discussed qualitatively 
in order to exhibit the emerging structure. Owing to the complexity of this procedure, only terms proportional to 
the third power of the electron-phonon coupling are found in Sect ion ITTT1 The resulting extended Frohlich Hamiltonian 
is transformed in the next section by a second Frohlich- type transformation, which produces 4-electron terms. These 
are discussed in detail in Section In particular, a reduction similar to the BCS one is performed, which yields 
an interaction of the form Vmt in Eq. The expression for gkk' is derived. Detailed analysis of this expression 
is performed in the next section, in particular, we showed, by applying some approximations, that this expression is 
negative i.e., Vmt is attractive. The original Frohlich's transformation and the resulting terms, up to 6th power of 
the coupling, are commented in Section IVHI The final section contains a discussion, summary and open questions. 



II. HIGHER ORDER TERMS OF THE FROHLICH'S TRANSFORMATION 



Following Frohlich « (see Appendix 1X1 for details), let us consider the electron-phonon Hamiltonian: 




& wPw) 



(4) 



where 




(5) 



and Oko- (&k) are fermion (boson) operators. The coupling D w will be assumed small and % = 1. 

Since the interaction is spin-independent, the spin index will be suppressed. Summation over electron momenta 
will include summation over spins. 
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Frohlich performed a unitary transformation of H e - p h in order to eliminate (as far as possible) the interaction term. 
The transformed Hamiltonian is 

H = e s * ff c _ ph e s = H c _ ph - [S, ff - P h] + \ [S, [S, H c _ ph }} + . . . , (6) 

where 

s = E s « = E - ^) = ~ s *> ( ? ) 



7q = E ^ <l) a k a k-q, (8) 
k 

and the unknown function (/>(k, q) : R 3 x R 3 — + C 1 is adjusted to achieve the cancellation. 

Subsequently, a term which is a combination of products, each with / fermion operators and b boson operators will 
be written as (f,b). Clearly, / will always be even. For example, Hq consists of terms (2,0) and (0,2). 

The rhs of Eq. © expresses in terms of commutators [(fx, &i), (/2, 62)]- O ne easily finds that 

61), (/a, 62)] = [fuHbxh + hhlbxM = [/i,/a]W>i + /i/a[&i,&a]. (9) 

The necessary commutators [fi,/?], [61,62] are given in Appendix IO 

According to Eq. ©, the transformation can be performed, given commutators of the form occurring in Eq. (JSJ 
with the first argument equal S. The latter is a (2, 1) expression, hence 

[S, (/, b)} = [(2, 1), (/, 6)] = (/, b+ 1) + (/ + 2, b - 1), (10) 

by virtue of Eqs. ©, GH- Clearly (/, b - 1) = for 6 = 0. 

Based on this ground, Frohlich obtained the transformed Hamiltonian (see Appendix |A"|) : 

rr V- 1 \ ^ £%(1 + A(k,w))(l- A(q,w)) . , 

H F = 2^ £ k«k - r 2^ _ c — {a^a^a^a^ + c.c.) . 

k kqw 

The second term represents an effective interaction between electrons dressed in the phonon field. If e q - w — e q +w w > 0, 
this interaction is attractive. 

Omission of higher terms in Eq. ||f)J results in violation of unitarity. The question thus arises whether partial 
inclusion of these terms (first of all those of 3rd order in the coupling) could improve agreement between theory and 
experiment. 

Following the general rule for the action of S in consecutive orders, expressed by Eq. I|10|) . one easily finds the form 
of subsequent terms: 

{£S} A ^ A {£S} A ^ A {£S} A ^ A {£S}- 

In each consecutive step one obtains terms proportional to the next power of the coupling. 

From the view-point of Fermi liquid theory, the terms representing effective inter-electron interactions are most 
interesting. According to Eq. (|llfl . new terms (6,0) (proportional to D^), describing 3-electron interactions and (8,0) 
(proportional to -D^,), representing 4-electron interactions, appear. Attractive interactions of this type leading to 
the formation of fermion triples and quadruples, could affect the behaviour of a superconductor. However, the total 
spin of an electron triple is nonzero, so such clusters are unstable, as they are not invariant under time inversion. So 
far, there has been no experimental evidence of such objects. 

Most electrons in a superconductor below T c are paired, so the 4-electron interaction between Cooper pairs can be 
expected to prevail. Furthermore, quadruples with a total spin equal zero and appropriate 1-electron momenta are 
stable under time inversion. On the other hand, under Frohlich's conditions for convergence of series 0, the effect 
of the terms (6, 0), (8, 0) is weaker. 

Evaluation of these terms is a complicated procedure. Before doing this, let us first examine the 3rd order corrections. 
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III. THIRD ORDER OF THE TRANSFORMATION 

Let us consider the effect of the first higher orders discarded by Frohlich, i.e., terms proportional to D^. Then 
the corrected Hamiltonian takes the form 

H'=H Q - ([S, H ] - H mt ) + Q [S, [S, H }] - [S, H M ]j - Q [S, [S, [S, H ]]] - \ [S, [S, H mt ]]^j +.... (12) 

The additional two terms in the last bracket are equal, explicitly, 

\ [S, [S, H int ]] - \ [S, [S, [S, Ha]]] = £ A kq 6 q n k7q 



qk 

^ { ^kqwk' fflk-w fl k-q a k' -q fl k' + Ckqwk' O^' ak'-qOk-qOk+w } + C.C. != H e , 



(13) 



qkwk' 

where 

2A kq = iDq<j)*(k, q) + iD^>(k + q, q) + ~ (e k - q - £ k + w q ) |0(k, q)| 2 - i (e k - e k+q + w q ) |c/>(k + q, q)| 2 + c.c, (14) 

2S kqwk , = iD^* (k', q)0* (k - q, w) - 0* (k', q)0* (k, w) + 0* (k, w)0(k, q) - 4>* (k - q, w)0(k - w, q) } 

+ -^*(k',q)|(/)(k, q)0*(k, w) (e k - e k _ q + £ k -_ q - £ k ') + </>(k - w, q)0*(k - q, w) (e k _ w _ q - e k _ w - e k /_ q + £ k ')j, 

(15) 

2C kq wk' = iD q {0(k', q)0*(k + w - q, w) - 0(k', q)0*(k + w, w) 
- 4>* (k + w - q, w)0* (k + w, q) + 0* (k + w, w)(f>* (k, q) j 
+ i</>(k',q){V(k + w,q)0*(k + w - q, w) (E k '-q - £k' - £k+w- q + £k+w) 
+ 0*(k, q)0*(k + w, w) (e k _ q - e k - e k /_ q + e k ')|, 



(16) 



It can be seen that a new phonon index w has appeared in Eq. (|13|) . It results from the commutator of S with 
the terms (4, 0), so the harmonic approximation has not been violated. 
Substitution of the expression (|13l) into H' yields 

H = H a + Hb + H c + Hd + H e , 

where 



^ab*b a + c.c. 



H b = ^6*6 q 7i k |i£) q (0(k + q,q) - 0(k,q))-i (e k - £ k+q + cj q ) |0(k+q, q)| 2 + i (e k _ q - e k + w q ) |0(k, q)| 2 }+c.c, 

qk 

(17) 

H c = ^2bq(iD cl + (£ k _ q -£ k + w q )0(k,q))a k a k _ q + c.c, (18) 
qk 

H d = ^2 f^q^lk-q) + ^(£k- q -£k + ^q)0(k,q)0*(k',q) j a k a k _ q a k ,_ q a k / +c.c.. (19) 

qkk' ^ ' 

Hb is linear in 6 q 6 q but differs from H a , H c is part of Hi nt which had been excluded from the transformation to avoid 
convergence problems, Hd is the 2-electron interaction obtained, whereas H e contains terms of the form (4, 1) which 
represent the obtained correction to Frohlich's Hamiltonian. 

The function 0(k, q) remains in the form (|A7|) . which guarantees minimization of the contribution H c . 
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IV. THE SECOND TRANSFORMATION 

The structure of H' bears similarities to that of H e ^ p h . Both H' and H e ^ p h contain terms describing the interaction 
of electrons with phonons: the counterpart of H int in H e - P h is H e in H' . To estimate the effect of H e , let us repeatedly 
apply Frohlich's method and perform a second unitary transformation adjusted to eliminate H e as far as possible. 
The form of H e suggests to take 

s' = J2 « 6 u - vub u + £X - Cu&u + CK - M , 

u 

where 

Vu = Yl m ' u Kaia m _ u a m , (20) 

lm 

C = ^2x*(lm,t,u)ai_ u a^ t a* m _ t a m , (21) 

lmt 

C = ^^*(l,m,t,u)a^ai_ t am-t a m+u. (22) 

lmt 

The explicit form of the functions ip, x, f will be found below. 

The interaction H e is 3rd order in D q , therefore bearing in mind Frohlich's approach, we shall restrict the expansion 
of exp[S"*]_ff' exp[S"] = H to terms which are 6th order in _D q . These terms are 

H = H'- [S', H'] + i[5', [S\ H a }] + .... (23) 
To evaluate the rhs, one needs the commutators 

[£", H a ] = ^2 K^*( 1 ' k , u ) ( £ k - £k-u - w u ) a!*aiak_ u a k 

ukl 

+ ^2 b* u X* (1) k, t, u) (e k - e k _ t + £i-t - ^l-u - <*> u ) a^_ u a^ t al_ t a k 

uklt 

+ X! ^^(Ik^u) (£k+u - £k-t + ei-t -£i - w u ) ajai_ t a k _ t a k+u + c.c, 

uklt 

[s\ H b ] = -J2 « + q + c*) (^ kq + ^ kq ) + c.c, 

qk 

where D kq is defined by Hb — ^ q o*6 q ?i k £> kq + c.c. in Eq. I|17|). 

#c] = - J] ^k.qafcflk-q « + + Cq) 
qk 

+ X^ k , q 6*6 q |^*(l,k,q)ni (n k _ q - n k ) + (^*(k,l,q) - ^*(k - q, 1, q))a k a k _ q ai*_ q ai} 

qkl 

+ X 6 q 6 q{(- S k,qX*( 1 'M 1 q) - _E k _ t ,qX*(l,k- q,t,q))a 1 *_ q ai_ t a k _t a k- q 

qklt 

+ (#k+ q -t, q X* (k + q, 1, t, q) - Sk,qX* ( k , 1, t, q)) a k a k _ t ai*- t ai 

+ (-E k+q+t , q ^*(I,k + t,t,q) - B kiq ^*(l,k + t - q, t, q))a 1 *ai_ t a k ak+t 

+ (-B k , q <^*(k + t,l,t,q) - E k+t ,q(p* (k + t - q, 1, t, q))a k+t a k _ q a 1 *_ t ai +q \ + c.c, 
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where i?k,q is defined by H c = X) q k ^q a k a k-q-^k,q + c.c. in Eq. JTSJ, 

[5', = ^ 6* (F k)k / )q + F^ ikiq ) |a k a k _ q [77* , a k ,_ q a k ,] + [7/* , a k a k _ q ] a k ,_ q a k , 

uqkk' 

+ a k ak-q [£*, a k '-q a k' ] + [C a k ak- q ] a k /- q a k ' + a k flk- q [C a ik'-q a k'] + [C, a k a k - q ] a k ,_ q a k ,| + c 
where F k , k ', q is defined by H d = J] qkk ' F k;k / jq a k a k _ q a£.,_ q ak/ + c.c. in Eq. O, 

[5', ffe]=-E A ^K b *{ ^7q] + [Cq, «k7q] + [Cq, ™k7q] } 

qk 

— ^ Sk,q,w,k'^w^w| [»7w) ak-w a k-q a k;'-q a k'] + [?w, ak-w a k-qak'~qOk'] + [Cw , Ofc- w a k-qO k ' -q a k'] | 
qwkk' 

— ^2 Ck,q,w,k'^w^w| [Vvfi ak' a k'-qak;-q a k+w] + [? w , Gfc' fl k' -qOk:-q a k+w] + [Cw , Ofc' fl k' -q a k;-q a k+w] J 
qwkk' 

— ^k,q{?7q"k7q + Cq«k7q + Cq"k7q} 
qk 

— ^ ] Sk,q,w,k' |^^w a k-w a k— qOk'-q fl k' + Cw a k-w fl k-qOk' — q a k' — w 
qwkk' 

— ^ C k ^ q , W! k'|f7wa k /a k /_ q a k _ q a k+w + ^ w a k 'a k '_ q a k _ q a k+w + Cwa k /a k /_ q a k _ q a k+w | + c.c, 

qwkk' 

where A kq , £?k, q ,w,k', Ck, q , w ,k' are given, respectively by Eqs. Q , i jTB|. 

[S", [S",.ff ]] = - ^ G k ,k',q^ q ^q|[?7q,ak'ak'ak_ q ak] + [£q, a k 'ak'a k _ q ak] + [C q , ak' a k'a k - q «k] } 

q kk' 

- X] ^k,k',t, q 6 q fc q | [r/q, ak'- q a k'-tak-t fl k] + [C q > a k'- q a k'-tak-t a k] + [C q : a k'- q a k'-ta k _ t a k ] I 

cjkk't 

- J! ^k,k',t, q & q & q | [?7 q , a k ,a k /_ t a k -t fl k+ q ] + [C qi a k' a k'-tak-t a k+ q ] + [Cq> a k' a k'-tak-t a k+q] } 

q kk't 

- ^2 Gk,k', q |?7 q ^k'a k _ q a k + C q "-k'a k _ q a k + C q ™k'a k - q ak j 

q kk' 

- ^2 ^k,k',t, q |f? q ak'- q a k'-tak-t a k + ^ q a k ,_ q a k '_ t a k -t a k + C q °k'- q a k'-tak-t a k j 

q kk't 

^2 ^k,k',t, q jf7 q ak' a k'-tak-t a k+ q + C q a k' a k'-tak-t°k+ q + C q a k' a k'-ta k _ t a k+q | + c.c, 



(24) 



q kk't 



where 



G k ,k', q ='0*(k / ,k,q)(£ k -£ k _q-a; q ) ) (25) 
#k,k',t, q = X*(k', k, t, q) (e k - e k -t + e k '-t - £k'- q - w q ) , (26) 

^k,k',t, q = <y3*(k',k,t,q) (e k + q - e k -t + £ k '-t - £ k ' - w q ) . (27) 

The commutators in the terms containing phonon operators have not been evaluated, as we are interested first of all 
in expressions containing exclusively electron operators. 

Given the transformed Hamiltonian H , we are now in position to minimize the effect of H e by imposing, similarly 
as Frohlich, the condition H e — [S', H a ] = 0. This leads to equations for tp, \, (p, which determine these functions 
uniquely, viz, 

r(k',k,q)= (28 ) 

£k — £k- q — w q 
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After substituting <j> and A given by Eqs. l|A7)l . fR). one obtains 

^(k-,k, q ) = 2 ^ (1 - A(k -' 1) U 1 - A ' k '"' _i-^+i.<jy (29) 

3 (e k -q -£k + w q ) V £ k'- q - £k' + w q e k ' - £ k '+ q + w q / 

Introduction of additional functions in order to preserve convergence is not necessary here, as A already guarantees 
this property. 

As for x an d f, additional functions preserving convergence are indispensable, viz, 

X*(k',k,t,q) = J?K±5± (i _ A(k,k',t,q)) , (30) 

£k - £k-t + £k'-t - £k'-q -w q \ / 

where 

A(k,k',t, q) = {^ !! j^-^-t + ^-t-^- q -o; q |<r q _ 

(^0, it |£k-£k-t +£k'-t -Ek'-q-W q | > T q 

Taking into account Eqs. (|A7I) . I|15|) . one obtains 

x . (k ; k , t , q) = 2^ ( (i-^*)) (i -*(*■«)) + 

6 (£k - £k-t + £k'-t - £k'-q - W q j V (£ k -t - £k + LUt) (£k'-q - £k' + w q ) 
(1 - A(k, t) ) (1 - A(k' - t, q)) (2 + ^ _ ^ t)) + (l-A(k-,t)) (l-A(k-,q)) (4 _ A(k) t)) 



(£k-t - £k + W t ) (£k'-t-q - £k'-t + Wq) (£k'-t - £k' + W t ) (£k'-q - £k' + Wq) 

(l-A(k'-q,t))(l-A(k'-t,q)) 



(£k'-q-t - £k'-q + Wt) (£k'-t-q - £k'-t + W 



(4-A(k,t)) 1-A(k,k',t,q) 



(32) 

^*(k', k, t, q) = (l - A(k, k', t, q)) , (33) 



Ck,t,q,k 

£k+q — £k-t + £k'-t — £k' — W 

where 



A(k,k',t,q) = J 1 ' if j £ k + q-£k- t +£k'-t-£k'-a;q|<f q; 

[0, if |£ k+q - £k-t + £k'-t - £k' - W q | > T q 

and substitution of </> and Ck,t, q ,k' from Eqs. (|A7(I . fTHjl yields 

^ ( k', k, t, q) = - ^ ( (l-A(k-,t))(l A(k + q-t,q)) + ^ + 

6 (£ k +q - £k-t + £k'-t - £k' - W q ) V (£k'-t - £k' + Wt) (£k-t - £k+q-t + W q ) 

(1 - A(k', t)) (1 - A(k + q, q)) + (l-A(k + q ,t))(l -A(k + q-t,q)) _ 



£k'-t — £k' + ^t) (£k — £k+q + ^q) ' (£k+q-t — £k+q + ^t) (£k-t — £k+q-t + ^q 

(l-A(k,t))(l-A(k + q,q)) 



(4-A(k',t)) 1-A(k,k',t,q) 

(£k-t - £k + W t ) (£k - £k+q + Wq) /V 

(35) 

Having established tp, \, ip, let us average H over the phonon vacuum: 

H av = g XI £knk + X ^k,k',qak a k-qak'-q a k' + Sk, q ak a k-q (^ q + £ q + Cq) 
k qkk' qk 

qkk 

I 

qkk 



^ Mk',q^*(k, q) - -Gk,k',q j |?/q«k'ak-q a k + £q«k'ak:-q a k + Cq"k' a]k-q a k } 
ikk' ^ ' 

X ( 5 k',t,q,k - 2^k,k',t,q J |??qak'-q a k'-ta k _ t ak + Cq a k'-q a k'-ta k _t a k + Cq a k'-q a k'-ta k -t a k| 
[kk't ^ - / 

^ f Ck,t,q,k' — --fk,k',t,q j |^ q a k 'ak'-tak-t a k+q + Cq a k' a k'-tflk-t a k+q + Cq a k;' a k'-tak-t a k+q| + C.C. 
ikk't ^ ~ ' 



(36) 
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H av contains the free-electron term, the 2-electron Frohlich interaction and terms representing 3-electron and 4- 
electron interactions. The 3-electron terms were generated by the second transformation of H c . The source of these 
terms is thus the non-transformed part of the original interaction, discarded by BCS. If higher-order corrections to 
BCS theory are of interest, all terms arising from that part, i.e., the 3-electron ones, can be therefore neglected. 
The same conclusion was drawn above on the grounds of unstability. 

V. 4-FERMION INTERACTIONS 

The Hamiltonian H av contains several terms representing 4-fermion interactions. Using Eqs. I|25[l. I|28|) . I|30|) . i|33[l . 
one finds 

^k', q 0*(k,q) - ^G k , k <, q = ^Gk.k'.q = ^*(k,q)A k ', q . 
Additionally, taking into account Eqs. i|2ti|) . 1|27|) . we get 

^k',t, q ,k — 2^ k k ' t i = 2^ k ' k '' t,q = o^ k ' ,t ' q ' t {} ~ ^0*-' ^)) ' 

Ck,t, q ,k' — 2^ k > k '>*>^ = 2^ k > k '> tq = 2^ k,t,q,k ' {} ~ ^O^k'jtjq)) • 
In terms of Gk,k',q, #k,k',t,q, -fk,k',t,q the 4-fermion interactions present in H av express as 

H\= ^G k , k ',qV'(l ) m,q)a^ 1 a m _ q nink'a k _ q ak + c.c, (37) 

qkk'lm 

Hf = ^2 ^ G k,k', q x(l, m, t, q)a^ 1 a m _ t a 1 *_ t ai_ q n k 'a k _ q ak + c.c, (38) 

qkk'lmt 

flj = 2^ k ' k ''^^' m ' *> ^ a m+q a m-toI-t a i n k'aJ_ q ak + c.c, (39) 

qkk'lmt 

i?4 = ^2 ^k,k',t,q^(l, m,q)aj n a m _ q nia k ,_ q a k '-ta k _t a k + c.c, (40) 

qkk'lmt 



#1 = X! 2 1 ^ k ' k '' t ^ x ^' m ' w ' q ^ a ™ am_wai *- wai_qak '-q ak ' _tak - tak + c ' c '' *- 41 ^ 

qkk'lmtw 

Hl= ^k,k%t, q ^(l, m, w, q)aj n+q a m _ w a 1 *_ w aia k ,_ q ak'-ta k _t a k + c.c., (42) 

qkk'lmtw 

H J = ^ ^k : k',t, q ^(l,m, q)aj n a m _ q nia k ,a k '-ta k -t a k+q + c.c, (43) 

qkk'lmt 

iff = 2^ klk '' t ^ qX '~ 1, m ' W ' ^) a m fl m-wflr-w a l-q a k' a k'-tak-t a k+q + C.C, (44) 

qkk'lmtw 

i?f = ^ -^k,k',t, q ^(l, m, w, q)a m+q a m _ w a 1 " 1 _ w aia^ak'-tak-t a k+q + c.c. (45) 

qkk'lmtw 
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TABLE I: Reductions of Hi 
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m — > m — 2k 


2 


k 


m 


-2m 


k — m 


k + m 


// /// / 

a — — cr = — cr = a 
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k + m 


k + m 


ii in i 
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m 
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a — — O" ,(7 = — O" 





TABLE II: Reductions of Hi 



Five of these interactions, viz, H\, iff, iJ|, H%, ifj, contain the operator rik, therefore they are not reducible to the 
4-fermion MT potential • This potential is particularly interesting not only for its possible relevance to the physics 
of superconductors, but also because the thermodynamics of the Hamiltonian Hbcs + W + Vmt is exactly solvable 
(Brankov et ali^i). 

Let us consider the 8-fold product of operators in fff with explicit spin indices: 

* * * * 

With respect to 1-fermion momenta these are nine possible reductions to the form ©, and each of them allows two 
or four possibilities related to spin reduction. The values assumed by each momentum index are collected in Table [I] 
In four cases q = and the coupling vanishes, since x(l, m > w > 0) = 0, for all 1, m, w. As a consequence, H\ assumes 
the reduced form 

red) — ^2 { 2 ^k,2k+m,k+m,2k+2m x( 2m + k, HI, k + HI, 2k + 2m) 

- #k,m,k+m,2m(x( 2m + k, m, k + m, 2m) + x(2m - k, m, m - k, 2m)) 

— x(k, m, k + m, 2k) (-ffk,2k+m,k+m,2k + #k,2k-m,k-m,2k) jam- a m+ a -m- a -m+ a -k+ a -k-ak+ak- + C.C.. 

(46) 

Terms with k = m have been excluded, similarly as in Vbcs- Their contribution is accounted for by a shift of 
1-fermion energies. The coefficients on the rhs result after performing summation over spins. 

This procedure has been also applied to H%, Hf H®. The corresponding values of 1-fermion momenta are given in 
Tables ITT1 IIIII IIVI The additional remark k — > k 2m (or similar), indicates that a translation of one momentum 
index is necessary after reduction. 
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a" = -a' = -a = a'" 




9 


-k 


m 





k m 


m k 


a" = -a', a'" = -a 





TABLE IV: Reductions of fff 



After reduction, U|, iff, iff take the following forms: 

H t(red) = { 2 #k,2k+m,k+m,2k+2m ¥>(-k, -2k HI, k HI, 2k + 2m) 

- #k.m,k+m,2m(<^(k, -m, k m, 2m) + (f(-k, -m, -m - k, 2m)) 

- ip(—\c, m - 2k, -k + m, 2k) (if k ,2k+m,k+m,2k + -ffk,2k-m,k-m,2k) |a^ n _a^ 1+ a , l m _a !, L m+ a_ k+ a_ k _a k+ ak- + c.c., 

(47) 

^4(red) = X! |2/-k-2m ,-m ,-k-m,2k+2m x(k + 2m, m, k + m, 2k + 2m) 

k^m 

- /k-2m,-m,k-m,2m (x(k + 2m, in. k + in, 2m) + x(-k + 2m, m, m - k, 2m)) 

- x(k, m, k + m, 2k) (l-k,m,-k+m,2k + -f-k,-m,-k-m,2k) |am- a m+ a -m- a -m+ a -k+a-k-a k +ak- + c.c, 

(48) 

H l{red) = ^2 { 2/ -k-2m,-m,-k-m,2k+2m ¥>(-k, -2k m, k m, 2k + 2m) 
k^m 

- 4-2m,-m,k-m,2m(^(k, -m, k - m, 2m) + (p(-k, -m, -k - m, 2m)) 

- ¥>(— k, -2k + m, -k + m, 2k)(l_k,m,-k+m,2k + ^-k,-m,-k-m,2k) } x a^ 1 _aJ n+ a!L m _al m+ a_k+a-k-ak+ak- + c.c. 

(49) 

Collecting all terms in Eqs. I|46|l - H49(l . one obtains a Vmt interaction as in Eq. © with 

3mk = -4w 2k +2mA km A^ k + 26 km + 26^ k , (50) 
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where 

Akm = x*(2k + m,k,k + m, 2k + 2m) + ip*(-m, -k - 2m, -k - m, 2k + 2m), (51) 
©km = w m (x*(m,k, k + m, 2m) + tp*(—m, k - 2m, k m, 2m)) 

x ^x(2m + k, m, k + m, 2m) + x(2m k, m, m — k, 2m) + </?(k, — m, k — m, 2m) + (p(— k, — m, — k — m, 2m)^ . 

(52) 

VI. 4-FERMION INTERACTION COUPLING 

The most important question about the 4-fermion interaction is positive- or negative- valuedness of g m k. In particu- 
lar, it would be of interest to determine the domain of g m k in momentum space where it is attractive. Unfortunately, 
the form of g m k given by Eq. (|50|) is extremely complicated, so this problem cannot be resolved in general. 

First, let us specify the quantities occurring in Eq. J50J. We assume that £k = ak 2 , Wk = bk, I\ = ck, Tk = dk and 
Tk = ek, where a, b, c, d, e are real, positive constants. 

Following Davydov^, we apply some approximations in order to estimate g m k- We assume that the most significant 
contribution to the 4-fermion interaction (analogously as in the BCS theory) comes from 1-electron momenta m, k 
which satisfy the condition m ~ k. 

The case of exact equality is, of course, excluded (in accord with the restriction on summation in Eqs. I|46|) - l|49|lh 
so we use gkk only as an abbreviation for g m k under our approximation m « k. 

Thus, the interaction coupling gkk is given by 

Skk = -4cj 4 kAkkA k k + 29 k k + 2e kk - (53) 

We take into account only vectors with almost compatible directions, otherwise gkk vanishes. 
We can thus rewrite Eq. (|53[) in the form: 

.g kk = -16^ k |A kk | 2 + 43fJ6 k k, (54) 

where the first term is always nonpositive and 5R means real part. The second term requires detailed analysis. Detailed 
calculations are performed in Appendix [B] We have found gkk in all cases for different values of the constants a, 
b, c, d, e. Most interesting is the case b > c, b > e and b > d, because the essential part of initial Hamiltonian is 
transformed under these conditions. Under a further approximation, we find the form of gkk for k ss kp, where kp is 
Fermi momentum. Additionally we put c = d = e and £k F S> Wk F (this is at least true for metals). Under these 
assumptions 



which implies 



QkFkF * 4< ' ~ 12< £k , 



ffk F k F ~ — § — < 0. 



Now we can compare the magnitudes of the coupling constants of 4-fermion interaction and 2-fermion interaction 
(under the same assumptions). It is shown in Appendix 1X1 that Gk F k F = --D^ F /wk F - Thus the interaction coupling 
for 4- and 2-fermion (BCS) interactions fulfils the relation (?k F k F = G^ FkF /uj^ F . As a consequence, for a strong 
pairing, there is a significant contribution from 4-fermion interactions. 

Moreover, the 4-fcrmion interaction coupling gkk is also negative in most cases without imposing the approxima- 
tions c = d = e and k w &f . This is shown in Appendix [Bj 



VII. 4-FERMION INTERACTIONS IN FROHLICH'S EXPANSION 



As demonstrated in Section llTTl 4-fermion interaction terms appear in sixth order of the expansion of Frohlich's orig- 
inal transformation. These terms are derived in Appendix El The various resulting 4-fermion terms considerably 
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outnumber those obtained by applying a second transformation in Section Again a reduction procedure to Vmt 
is possible for all 4-fermion expressions l|C5(l . which do not contain the particle number operator and possess three 
phonon indices. However, the resulting couplings are also complicated functions and, therefore, will not be examined 
in detail. 

The disadvantage of this method is appearance of 3-fermion interactions, which in the double transformation method 
gave a small contribution. Here additional arguments must be used to discard these terms. 

VIII. CONCLUDING REMARKS 

We have extended Frohlich's transformation of i/ e -ph to higher order terms. This has been done by performing 
a second transform of the first terms discarded by Frohlich. The resulting interactions are of 3- and 4-fermion type. 
The 3-fermion terms can be expected to be inessential in superconductivity because of their instability under time 
inversion. The 4-fermion terms are, in general, reducible to Vmt, a BCS-type 4-fcrmion interaction. The resulting 
4-fermion coupling is extremely complicated, but under reasonable approximations it is negative- valued. Moreover, 
for 1-fermion momenta in the neighbourhood of Fermi momentum, it has the simple form — D^ f / uj^ f = G^. FkF /oj^ F , 
where Gk F kF is the 2-fermion coupling. This fact implies that 4-fermion interactions are significant for systems with 
strong pairing and allows to estimate (relative to 2-fermion coupling and phonon energy at Fermi momentum) their 
magnitude. 

As in BCS theory, where the relation between the gap parameter A and coupling constant G of the BCS interaction 
allows to estimate the magnitude of G, it can be expected that detailed thermodynamics of Hmt will provide 
a relation between |<? m k| and other parameters of the theory, thereby allowing to estimate the magnitude of g m k- 
Another emerging question is Cooper's problem 2 ^ for a bound quadruple in the presence of Vmt or Vbcs + Vmt- 
Kamei and Miyake deal with this question in a recent work* 2 *. 

The double transformation has unveiled the structure of 3-, 4-, 5-fermion interactions. Since they are proportional 
to 4th, 6th, 8th power of D w , they are relatively weak, so inclusion of these terms, apart from 4-fermion ones, appears 
unjustified at present, although Schneider et alii^ suggest existence of both quadruples and sextets in some HTSC. 

The higher-order expansion terms of the transformed 7? c _ p h, including quadruple, sextet etc. interactions can be 
expected to reveal themselves in materials with extremely strong pairing correlations between spin 1/2 fermions, since 
the presence of strongly correlated pairs implies at least some kind of weak interaction between them. 

On the other hand, our extension of Frohlich's transformation shows that if phonon mediation exists in a supercon- 
ductor, the 4-fermion and 6-fermion interactions are always present as supplementary to the BCS one. Unfortunately, 
the Hamiltonian -ffecs + Vmt, although exactly solvable, leads to intricate mean-field equations. 

Our results obtained can be generalized in many respects. First of all, the process of averaging over phonon vacuum 
could be replaced by averaging over the phonon equilibrium state, which could be justified at higher temperatures. 
This would already lead to an additional 1-fermion term in Frohlich's Hamiltonian and modification of 2-fermion 
and 3-fermion terms in our method. Another extension would result by going beyond the harmonic approximation 
and including all products of phonon operators. Method of Bogoliubov 25 could be applied to "dangerous terms" 
(divergent), omitted in Frohlich's method. These questions will be dealt with in further investigations. 

APPENDIX A: THE FROHLICH'S TRANSFORMATION 

Let us recall Frohlich's method^, using in some details the more elegant approach due to Davydov 22 . 
Consider the electron-phonon Hamiltonian: 

iJ _ph = H + Hint = ekOko-Okcr + ^ ^w&w & w + i ^ D ™ ( fo w/Ow _ kw/°w) , 

ka- w w 

where 

Pw = ^ ] a k— wcr a k(T, 

k(T 

and Ok CT (bit) are fermion (boson) operators. The coupling £> w will be assumed small and h = 1. 

The form of the interaction term of H c ^ pn arises under a number of assumptions: the ions of the lattice move 
collectively, the coupling depends only on w and electrons interact only with longitudinal phonons for which oj w = ws, 
s denoting the velocity of sound. Our interest is focused on the behaviour of electrons, therefore variations of 
the phonon spectrum will be accounted for only through s. 
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Frohlich performed a unitary transformation of fl e - p h in order to eliminate (as far as possible) the interaction term. 
The transformed Hamiltonian is 

H = e s *H c - ph c s = fl c _ ph - [S, fl e - ph ] + 5 [5, [S, fl e _ ph ]] + ... , (Al) 

where 

s = Y 5 q = (^q & q ~ Tq & q) = ^q = Y ^( k > <l) a k a k-q, 

q q k 

and the unknown function <p(k, q) : K 3 x R 3 — > C 1 is adjusted to achieve the cancellation. 
Collecting terms of the same order in the coupling Z? w , one obtains Frohlich's expansion: 

fl = fl - ([S, H ] - flint) + [5, flo]] - [5, flint]) + ■ ■ ■ ■ (A2) 

Subsequently, a term which is a combination of products, each with / fermion operators and b boson operators will 
be written as (/, b). Clearly, / will always be even. For example, H consists of terms (2, 0) and (0, 2). 
The rhs of Eq. I|A2(1 expresses in terms of commutators [(/i, b\), (fa, 62)]. One easily finds that 

[(fi,h),(f 2 ,b 2 )} = [hJ2]bib 2 + hh[b 1 ,b 2 ] = [/i,/ 2 ]Mi + /i/ 2 h,6 2 ]. (A3) 

The necessary commutators [fi,f 2 ], [^1,^2] are given in Appendix IO 

According to Eq. (|A1|) , the transformation can be performed, given commutators of the form occurring in Eq. (|A3|) 
with the first argument equal S. The latter is a (2, 1) expression, hence 

[S, (/, b)} = [(2, 1), (/, b)} = (f, b+ 1) + (/ + 2, b - 1), 

by virtue of Eqs. (|Alj) . itTTTTl . Clearly (/, b - 1) = for 6 = 0. 

One finds 

[S, H ] = [(2, 1), (2, 0)] + [(2, 1), (0, 2)] = (2, 1), 

so the term arising from Hq in first order, has the same form as the one arising from fli nt in zeroth order. These 
terms are used to eliminate the interaction in fl e _ pn . 
Similarly, for fl; n t 

[S, fl int ] = [(2,1), (2,1)] = (2, 2) + (4,0). 

Frohlich additionally assumed that introduction of experimentally measured velocity of sound s allows to discard all 
terms containing fr 2 ,, (6^) 2 , 6 W 6 V , b^,b v ifw^v etc. In other words, all (/, b) terms for b > 3 and part of (/, 2) terms 
are neglected. As a consequence, the expansion IjAljl up to 2nd order reduces to^: 

[S, fl ] = ^2 bq (e k - e k _ q - w q ) (j>(k, q)a k a k _ q + c.c, (A4) 

kq 

[S, flmt] = Y ^qTqPq + iDqb^b^k, q) (n k - n k - q ) + c.c, (A5) 

q kq 

[5, [5,fl ]] = ^ (e k - q - e k + w q ) 6 q 6 q |0(k, q)| 2 (n k - n k _ q ) + (e k _ q - e k + w q ) 0(k, q)a k a k _ q 7* + c.c. (A6) 

kq kq 

is now adjusted so as to minimize the contribution of terms (2, 1). This is achieved by the choice 2 

<Hk,q) = ~ lD \ (1-A(k,q)), (A7) 

£k-q — £k + W q 



where 



A(k,q)= J< * |^q-^ + -q|<rq. (A8) 

10, if |e k _ q - e k + cj q | > T q 
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A(k, q) is introduced to avoid divergence of the series IjAlfl and T q is positive energy choosing for convergence. 
Eqs. mj|-|(lZJ| yield 



H = ^£ k ak a k + ^w&w^w + J^&w & w (n-k - «k-w) ( 2 (^k-w - £k + w w ) |0(k, w)| 2 - iD w <p(k, w) + c.c 

k w wk ^ 

, ,Vn ft . i * * u , , 1 ^ (1 + A(k, w)) (1 - A(q, w)) , . 
+ 1 2^ (o w a k a k _ w - 6 w a k _ w a k J A(k, w) - - ^ — (a k «k-wa q - w flq + c.c 



2 ^ e c 

kw kqw 



(A9) 



Discarding the non-transformed part of the initial interaction (4th sum on the rhs) and taking the average in phonon 
vacuum, one obtains the Frohlich's Hamiltonian: 

„ V- 1 ^ £%(1 + A(k,w))(l-A(q,w)) , . . > 

H F = £ k^k - ^ 2^ Z Z Kflk-wa q _ w a q + c.c.) . (A10) 

^ . £q— w T" 

k kqw 

The second term represents an effective interaction between electrons dressed in the phonon field. If e q _ w — e q +w w > 0, 
this interaction is attractive. 



APPENDIX B: CALCULATION OF g kk 

Our objective is to find the explicit form of the following expression: 

Skk = -16^ k |A kk | 2 + 43fJ6 k k. (Bl) 

Taking into account Eq. I)52[l. we have 

e kk = ^ k (x*(k, k, 2k, 2k) + ip*(-k, -k, 0, 2k)) 

x (x(3k, k, 2k, 2k) + x (k, k, 0, 2k) + <p(k, -k, 0, 2k) + p(-k, -k, -2k, 2k)) (B2) 
= io kX * (k, k, 2k, 2k) ( X (3k, k, 2k, 2k) + v>(-k, -k, -2k, 2k)) , 

because k, — k, 0, 2k) = x(k, k, 0, 2k) = ^(k, — k, 0, 2k) = 0, as can be seen from the explicit form e.g., of the first 
of these functions (see Eq. 133(1 '): 



<p* (-k, -k, 0, 2k) = fj^ (l - A(-k, -k, 0, 2k)) (1 - A(k, 2k)) 



6w k ^o 

x {(1 - A(-k,0)) (-4 + A(k,0)) + (1 - A(k,0)) (4- A(-k,0))} = 0, 



because A(-k, 0) = A(k, 0) = 1, for all k (see Eq. (A8j). Similarly for x(k, k, 0, 2k) and tp(k, -k, 0, 2k). 
Similarly, from Eq. 1(51(1 . we get 

A kk = X* (3k, k, 2k, 4k) + ip* (-k, -3k, -2k, 4k) . (B3) 
On the grounds of Eq. ((30() , we obtain 



- (1 - A( ^:^ k ' 2k)) + *™ - - - W 



(B4) 



The functions A and A are equal or 1, depending on their argument (see Eqs. l[3*T}l. l(A8(l). so we obtain several 
conditions. The most important one is 

A(k,k,2k,2k) = w k > fk & > c, (B5) 
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and x*(k, k, 2k, 2k) = otherwise. Other conditions are not so strong, because they do not destroy all of parts of 
Eq. fBlfl . We have 

A(k, 2k) = ^=> uj k > T k b > d, (B6) 



A(-k, 2k) = ^ |4e k + oj k \ > T k <=> fc > ^— -. (B7) 

4a 



Taking into account Eqs. (|B5I) . (|B6I) and (|B7|) we get 

• If 6 < c, then x* (k, k, 2k, 2k) = 0. 

• If b > c and b > d, then 

3w k (w k + 4e k ) 

• If b > c, b < d and fc G (0, ^] , then x*(k, k, 2k, 2k) = 0. 

• If b > c, b < d and k G oo), then 



X *(k,k,2k,2k)= ^ *^ 



8w k (w k + 4e k ) 

Proceeding similarly, we get the form of other functions occurring in Eqs. 1B2() and i|B3l) . For x(3k, k, 2k, 2k), we 
have: 

• If b < c, then x(3k, k, 2k, 2k) = 0. 

• If b > c, b > d and k G (0,^) U (^,oo), then 

X(3k, k , 2 .c,2 1 c) = ^-' 5 "--y-) . 

3w£ (wk - 4e k ) 

• If 6 > c, 6 > d and jfe G [^, Mi] , then 



X (3k,k,2k,2k)- "" k 



n/2£ k 
4w, 3 



If 6 > c, 6 < d and fc G [0, ^] , then x(3k, k, 2k, 2k) = 0. 
If b > c, 6 < d and jfc G (^, oo) , then 

X(3k,k,2k,2k) = 



8w k (u k - 4e k ) 2 



For ^(-k,-k,-2k, 2k): 

• If b < e, then <p(-k, -k, -2k, 2k) = 0. 

• If b > e, & > d and fc G (0, ^) U (^,oo), then 

p(-k, -k, -2k, 2k) 

• If b > e, b > d and k G [^, ^] , then 



iA/2j3 £k (^ k + i2e k ) 

3^K-i64) • 



y>(-k,-k,-2k,2k) = - 

'■' v k 



^3 



k 
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If b > e, b < d and fc £ [0, , then p(-k, -k, -2k, 2k) = 0. 
If b > e, b < d and fc £ [^,oo], then 

<^(-k, -k, -2k, 2k) 



8w k (ojI - 16e£) ' 

For x*(3k, k, 2k, 4k): 

• If b < c, then X * (3k, k, 2k, 4k) = 0. 

• If b > c, b < d/3 and fc £ (0, 4j^] , then x* (3k, k, 2k, 4k) = 0. 

• If b > c, b < d/3 and k £ (4=^, 4^), then 

X* (3k, k, 2k, 4k) = — - — := X (1) • 

V ' 16c^ k (wk + 2e k ) K + 4e k ) 

• If b > c, 6 < d/3 and fc e (4^,oo), then 

^•' 3k ' k - 2k ' 4k »' 4 W -^K-44) :=x "' 

• If b > c, d/3 < b < d and fc £ (0, 4=^] , then X * (3k, k, 2k, 4k) = 0. 

• If 6>c, d/3<6<dandfc€ (4=^,4^], then x* (3k, k, 2k, 4k) 

• If b > c, d/3 < b < d and fc e (4^,oo), then x* (3k, k, 2k, 4k) = x (2) ■ 

• If 6 > c, d < b < 3d and fc e (0, ^) U (4±4, oo) , then 



y*(3k k 2k 4k) - {lUjl ~ 164) - ^ 

V C ' 4k, -6^(^-164)(^-4^ ) -^ ' 



If 6 > c, d < 6 < 3d and fc e [^, 4^), then 



X*(3k,k,2k,4k) = -^(-y + 404 + 22. k£k ) 



If 6 > c, d < 6 < 3d and fc e [4^, 4±^] , then 

x* (3k, k, 2k, 4k) = - 2 r k r7w^ K ' „ ^ : = X 



i£> k (3o> k + 4e k ) -= v(B) 
24w£ ( Wk + 2e k ) (w k + 4e k ) 



. If 6 > c, 6 > 3d and fc e (0, ^4) U (^, 4^4) U (^,oo), then X * (3k, k, 2k, 4k) - X (3) . 
. If6>c, 6>3dandfce then X *(3k, k, 2k, 4k) - X (4) ■ 

• If6>c, 6>3dandfce [4^, 4^] , then x*(3k, k, 2k, 4k) = X (5) • 
For p*(-k,-3k,-2k,4k): 

• If b < e, then (-k, -3k, -2k, 4k) = 0. 

• If 6 > e, 6 < d/3 and fc e (0, 4j^] , then ip* (-k, -3k, -2k, 4k) = 0. 
If b > e, b < d/3 and fc e (4^, 4£?] , then 



p*(-k,-3k,-2k,4k) = — — , ^-7 — :=<^ (1) . 

y V ^ 16c^ k K + 2e k ) K - 4e k ) v 



If b > e, b < d/3 and k £ (4^,oo), then 



If 6 > e, d/3 <b < d and k G (0, ^] , then <p*(-k, -3k, -2k, 4k) = 0. 
If b > e, d/3 < b < d and k £ (^, 4j^] , then <^*(-k, -3k, -2k, 4k) = <pW. 
lib> e, d/3 <b < d and fc £ (^r.oo), then <^*(-k, -3k, -2k, 4k) = <^ 2 ). 
If 6 > e, d < b < 3d and k G (0, ^) U oo), then 



^*(-k,-3k,-2k,4k) = 

If 6 > e, d < b < 3d and k G ^) , then 

^*(-k,-3k,-2k,4k) = 

If b > e, d < b < 3d and k G ^] , then 



KK-44)K-i64) 



— «( 3 ) 



-1D| (3^ + 14e k g; k + 40e k ) _ (4) 
48w£ (yg - 44) ( Wk + 4e k ) ! ~ ^ 



<^*(-k,-3k,-2k,4k) = '' Dk 



16wj| ( Wk + 2e k ) 



If b > e, d < b < 3d and k G then 

<^*(-k,-3k, -2k, 4k) 



( 3cJk _ 4£k ) 
24w£ ( Wk + 2e k ) {u k - 4e k ) 



-,(5) 



• If b > e, 3d < b and jfe G (0, ^) U ^) U (^,oo), then <^*(-k, -3k, -2k, 4k) 

• If b > e, 3d < b and k £ , ^] , then p*(-k, -3k, -2k, 4k) = <p( 4 ) . 

• If b > e, 3d < b and jfc G , ^] , then p*(-k, -3k, -2k, 4k) = <p( 5 ) . 

On the basis of this results and Eq. I|B1|) . we get the 4-fermion interaction coupling g kk : 

• If b < c and b < e, then g kk = 0. 

• If b < c, b > e, b < d/3 and fc G (0, 4=^] , then g kk = 0. 

• If b < c, b > e, b < d/3 and k £ (4j^, 4±^] , then 

— D 6 

■9kk = ; -j- -2 := ffkk < °- 

16w k (w k + 2e k ) (w k -4e k ) 



If 6 < c, 6 > e, 6 < d/3 and k £ (4±4\oo), then 



-£> k u k e£ (2 ) 

3kk = n ^2^/2 1fi 2^ : = ffkk < - 

( w k- 4e k) K- lfe k ) 



If 6 < c, b > e, d/3 < b < d and fc G (0, 4±&] , then g kk = 0. 
If 6 < c, 6 > e, d/3 < b < d and k £ ^] , then g kk = ,g k 

If b < c, b > e, d/3 < b < d and k £ (4j^, oo) , then 5kk = 5kk } 
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If b < c, b > e, d < b < 3d and k E (0, ^) U oo) , then 



4g|4K + i64f _ (3) 



5kk '^l-Aelf^l-l^lf 



If b < c, b > e, d < b < 3d and k e then 



-g£ (3^ + 14 £k ^ k + 404) 2 (4) 

ffkk = i/i/i 3/ 2 7^7 ~T^2 : = -9kk < o. 

144w k (w k - 4e k ) (w k + 4e k ) 



If 6 < c, 6 > e, d < 6 < 3d and k e , then 



3kk = — — - -2 < 0. 

16w k (w k + 2e k ) 

If 6 < c, 6 > e, d < b < 3d and k € then 

-Dj (3^ k - 4e k ) 2 (5) 
36w k (w k + 2e k ) (w k -4e k ) 

If b < c, b > e, 3d < b and k e (0, ^) U ^) U (^,oo), then ,g kk = 

If b < c, b > e, 3d < b and k E , then 5kk = ff <£. 

If b < c, 6 > e, 3d < 6 and k E ^] , then 5kk = 5^. 

If 6 > c, 6 < e, 6 < d/3 and k E (0, ^] , then ffkk = 0. 
If b > c, b < e, b < d/3 and k E ^] , then 



5kk = < °- 

8co k (w£ - i 6e 2) 2 



If b > c, 6 < e, 6 < d/3 and fc e (4=^, ^] , then 

-3Dj {ujj + 8e k ) (6) 

5kk = 77 } ,2/ 2 1fi 2,2 : = -9kk < 0. 

16w k (w k + 2e k ) (w k - 16e k ) 

If 6 > c, 6 < e, b < d/3 and fc E (4^, 00), then 



_ -g° K + 64^4 + 164) ._ m r „ 

■ 9kk ~~ o ^2 it 2,2 / 2 ^ 2,2 ' _ ff kk ^ U ' 

8w k K - 16£ k) K- 4 4) 



If 6 > c, 6 < e, d/3 < b < d and k E (0, 4=^] , then ,g kk = 0. 
If b > c, b < e, d/3 < b < d and k E (4=^, 4±&] , then 



.9kk = ; %- -3 := < 0. 

16w k (w k + 2e k ) (w k + 4e k ) 



If b > c, b < e, d/3 < b < d and k e (4$, , then .g kk = g™. 
If b > c, b < e, d/3 < b < d and k e (^,00), then 5kk = 5 ^ k } . 
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If b > c, b < e, d < b < 3d and k G (0, ^) U (^,00) , then 

4g|f| + 464^4 - 2848^4 + 46084) (9) 

■ 9kk " 9- k K-44) 2 K-i64) 2 • 5kk ' 

(9) 



We see that g^ < iff 



luI + 464c44 - 2848^44 + 46084 > 0. 
'.2 



Taking into account w k = bk and £k = ak , we obtain a 3rd order algebraic inequality 

4608Z 3 - 2848:cZ 2 + 464a; 2 / - x 3 > 0, 
for I = a 2 k 2 , x = b 2 . This can be easily solved^ — there exist only one real root ~ 0,002x, so 

fkk < iff k > k o } ~ 0.047&/a < b/Aa. 



If b > c, b < e, d < b < 3d and k € then 

„ kk = Dl 

144^ (^_44) 2 ( Wk + 4 £k ) 2 

x (9c4 + 348w k e k + 1396c44 - 2080c44 - 7616t44 + 7680w k e k + 184324) := 5kk 
It turns out that <? kk < for all k > 0, because (/ = ak) the equation 

18432/ 6 + 7680W 5 - 76166 2 / 4 - 20806 3 / 3 + 13966 4 / 2 + 3486 5 / + 96 6 = 
has only two real roots and both are negative: l 1 ^ ^ » — 0.208& and l[ 10 ^ — — 0.029&. 



If b > c, b < e, d < b < 3d and k G ^] , then 



£>£ (9w 4 + 144a£e k + 784^4 + 1632^4 + 11524) 

9kk — : t 9 < 0. 

36c4 (wk + 2e k ) K + 4e k ) 2 



If b > c, b < e, d < b < 3d and k G ^] , then 

= D| 

3kk 36w k ( Wk + 2e k ) 2 (a£ _ 164) 2 



D 6 

■9kk = ; ~w- 3 < 0. 

4 K + 4e k ) 3 (wk - 4e k ) 2 



(11) 



x (9(4 - 48w k £k - 832(44 - 2944(44 + 1664(44 + 18432w k 4 + 184324) := ff kk 
Skk 5 < iff fc e (0, 0.0736/a) U (0.4126/a, 00). 
. If b > c, b < e, b > 3d and k G (0, U U (^,oo), then 5kk = g®. 

. If b > c, b < e, b > 3d and k € , then 5kk = t^ '. 

. If 6 > c, 6 < e, b > 3d and fc G ^] , then 5kk = g<£> . 

• If b > c, b > e, b < d/3 and fc G (0, , then ffkk = 0. 

• If b > c, b > e, b < d/3 and fc G (^, then 



If b > c, b > e, b < d/3 and fc e [4=^, ^] , then 



= Dj K + 4 ^£k + 2£ k ) = (12) 

' 9kk 2 (o£ - 164 ) 2 + 2e k ) 2 (a* + 4 £k ) ' Skk 



If b > c, b > e, b < d/3 and fc e (^,00), then 



= pj H 4 ^4 + iw| + 164) = (13) 

5kk 4 K - 164) 2 K - 4e k ) 2 K + 4 £k ) • 5kk • 



9vk < for all fc > 0. 




If b > c, 


b > e, 


d/3 < 6 < d and fc e (0, ^] , then ,g kk = 


0. 


If b > c, 


6 > e, 


d/3 < b < d and fc G ^] , then 5kk 


(8) 

= Skk- 


If 6 > c, 


6 > e, 


d/3 < b < d and fc e 4£r] , then 3kk 


= & a) 


If 6 > c, 


6 > e, 


d/3 < b < d and fc e (^,00), then ,g kk = 




If 6 > c, 


6 > e, 


d < b < 3d and fc e (0, ^) , then 





32^4 



5kk 9c k - 164) 2 {u* - Aelf (u, k + 4 £k ) 

x (-3*4 - 20w k e k + 8444 + 40044 - 568w k4 - 2400^44 + 1152w k 4 + 46084) := 

fkk ) < iff fc G (0.2686/a, 00). 

If b > c, b > e, d < b < 3d and fc e then 

16£> k e k (60^ k + 225^ k e k - 120^ k e k - 752(44 + 960^ k e k + 23044) (15) 
9w k (<4_44) K + 4£ k ) 

9kk < for a11 k > °- 

If b > c, 6 > e, d < 6 < 3d and fc e ^] , then 

D k (81<4 + 1320w k £ k + 6544^4 + 13056w k £ k + 92164) 

.9kk — z , 5 < 0. 

9loI K + 4£ k Hw k + 2£ k ) 2 
If b > c, b > e, d < b < 3d and fc e ^] , then 

5kk 



9u£ (<4 - 164) 2 (w k + 4£ k ) K + 2£ k ) 2 
x (9c4 + 36w k £ k - 336w k £ k - 2560w k £ k - 3264t44 14 592c44 + 46080w k 4 + 36864£ k ) 

9uu < iff fc G (0, 0.1426/o) U (0.3456/a, 00). 

If b > c, b > e, d < b < 3d and fc e (^,co) , then g kk = . 

If b > c, b > e, b > 3d and fc e (0, ^) U ^) U (^,00), then 5kk = . 

If b > c, b > e, b > 3d and fc e ^] , then 5kk = g<£> . 

If b > c, b > e, b > 3d and fc e ^] , then 5kk = . 
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APPENDIX C: HIGHER ORDER TERMS OF FROHLICH'S TRANSFORMATION 

Evaluation of successive commutators of the expansion IjAlfl is a time-consuming task. Below the final results are 
given for the relevant orders. 

Since S = B — B* in Eq. 0, all expansion terms have the form g + g*. We focus on the first summand. 
The following formulae for boson operators 

[6i,6 2 6 3 ] = [6i,6 2 ] 63 + 6 2 [61,63] , 
[6162, 6364] = 6163 [6 2 , 64] + 61 [6 2 , 63] 64 + 63 [61, 64] 63 + [61, 63] 6462, 

and fermion operators 

[fi, hfs] = {fi, f2\ h — $2 i/i) h} , 

[/1/2, 73/4] = {/1, /3} fif-i - /3 fi} h + /1 {/2, h} h - hh {/2, fi} , 
will be used. In particular, for fermion creation and annihilation operators, 

[a* k ai,a*a r ] = Si q a* k a r - 5 kr a*ai. (CI) 

1. Fourth order 

To calculate 4th order terms, let us rewrite the 3rd order expression (113(1 as a sum of = i [S, [S, H int ]\ 

and^ 3) = -i[5, [S, [S,H ]}]: 

rr(3) \ A Aint,0r* * 1 \ A 7* f T-yint.Q * * . /~iint.Q * * j 

H U,0 = ^kq 0q"k7q+ 2^ °w ( 5 kqwk' a k-w«k-qOk'-q a k' + <- k qwk' °k' ^k' - q a k -q a k+w j + C.< 

qk qkwk' 



with 



A™ = ~i£> q (<£* (k, q) + 0(k + q, q)) + c.c, 



^kq = g (£k-q-£k + ^q) l<?Kk,q)| 2 - ^ (e k - £ k +q + Wq) |0(k + q, q)| 2 +c.c, 



^kqwk' = ^A, (0* (k', q)^* (k - q, w) - 0* (k', q)0* (k, w) + cf* (k, w)0(k, q) - 0* (k - q, w)0(k - w, q) 

-Bkqwk' = jU* ( k ', q) {</>(k, q)<?(>* (k, w) (e k - e k _ q + e k <_ q - e k >) 

+ (j)(k - w, q)0* (k - q, w) (s k _ w _ q - e k _ w - e k ,_ q + £ k - ) } , 

C k qwk' = \iD^ (^(k', q)0* (k + w - q, w) - 0(k', q)</>* (k + w, w) 
- 0*(k + w - q, w)0* (k + w, q) + 0*(k + w, w)0*(k, q) ) , 



C k qw k ' = ^(k',q)|0*(k + w,q)0*(k + w-q,w) (e k '- q - e k ' - e k +w- q + e k +w) 
+ </>* (k, q)0* (k + w, w) (e k _ q - e k - e k < q + s k > ) j . 



The 4th order term in Eq. JBJ equals 



H {i) = ~[S, [S, [S, H int ]}} + i [5, [S, [S, [S, H }}}} = -± 



Q rr(3) 



5,W .- ii int + H Q . 



Explicitly, 



rr(4) \ " .(4) „, \ ,-,(4) * * \ s*i(4) * * 

ff int,0 = 2^ Ant,o7q"k7q + 2^ B int : o7wa k -wak- q a k /_ q a k / + ^ ^int.olwak'ak'-qak-qflk+w 



kq kqk'w 



kqk'w 



X] Ant,0 & q 6 q"krik' + J! ^intVq 5 q a k a k-q7 q + F iit,0 & w fe w ffl k-w a k-q o k'-q ffl k' 



where 



G int,Q & w & wak'+w a k'-qak-q a k+w + ^ Ant.O^w^wakak-qflk'-q ^ + C.C. 
kqwk' kqwk' 



4 (4) — -4 



1 



int 



.4 



(4) 








B 



(4) 



is, 



int 



Mnt g^kqi ^0 ^^kqj '"int g^kqwk'i 

Other coefficients in Eq. (|C2|) arise according to the same scheme, viz, 



B, 



(4) 



-B 



o 

kqwk' • 



C$ = ^kqwk', ^ = i<(|0(k',q)| 2 - |0(k' + q,q) 



E$ = I (4£ q , q - O 0(k, q), F$ = l (B^ wM ,_„4>(V q, w) - 5L nt q, w ,k«/>(k', w)) , 



G$ = 5 (C^, w , k ^(k' + w, w) - C k n qiW!k , +w <Xk' + w - q, w) ) , 



4ft } = I (^ q) w,k^(k, w) - B k tw,q,w,k^(k + w - q, w) + a k f_ W)q)Wik ^(k' - q, w) - C^ q)W)k ^»(k' + w 



and similarly for Cq 4 \ D q 4 \ etc. 



2. Fifth order 



Proceeding analogously as with 4th order, one obtains 



[S, [S, [S, [S, [S,H } 



1 r 
4 



l r 

5 



5, H 4 ^ 



rr(5) , rr(5) 
^int + H 
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H b2,o = H ^mio 6 w7q«kak'-w a k'+ q + S £o 6 w7 q a k ak-w7q + Y c 'tot,o i C a k'-wOk'-qnk7q 

qkwk' qkw qkwk' 

-^nt,0 & u7wa k -w a k- q ak'-q a k'+u + ^ Ant,0 6 u7wa k -w-u a k- q ak'-q a k' 
qkwk' ii qkwk' u 

+ Y F int!o fo u a r-u a i-wa k -w a k-qak'- q ak' + Y G int,o fo u7wak- w ak-qa k '- q ak'-u 

qkwk'ul qkwk' u 

J_ \ ^ 7-(5) J * * | \ ^ x(5) 7 * * * 

T 2^ - / int,0° u ^ wa k-w a k-q-ua k /_ q ak' "h ^ j int,0° ua l+w a l-uttk_ w ak- q a k /_ q ak' 
qkwk'u qkwk'ul 

+ ^2 ^int ! o fe u7wa k /a k /_ q a k _ q _ u ak+w + ^inlo fe u7w«k'-u a k'-qak- q ak+w 

qkwk'u qkwk'u 

+ H M intVu a r-u a i-q a k ' a k'- q a k _ q a k '+w + 51 ^i ( nt,o 6 u7wa k 'a k '- q a k _ q a k +w-u 

qkwk'ul qkwk'u 

+ H °S,o 6 u7wa k , +u a k '- q ak_ q ak+w + Y P intV ua i+w a i-u<' a k'-q<_ q ak+w + Y i? int,o 6 q n k n k'7q 

qkwk'u qkwk'ul qkk' 

+ H S 'i ( nt,0 O q a k a k-q7q7q + ^ 5 'int,0 & q7qakak-q7q + ^ ^So^k- w^-q^ _ q a k ' - w7w 
qk qk qkwk' 

Y a int,o' , w fl k' + w a k'-q4-q ft k+w7w+ X! ^int!o 5 wa k ak- q a k '- q a k'7w + C -C-, 



qkwk' qkwk 

where 



42,0 = - { (42,0 + 4£o) >< (V (k', w)0* (k' + q, q) 0* (k' + q, w)0* (k' + q - w, q)) , 

o(5) _ _ f 1 H / ,(4) .(4). _ ,(4)k^k-w _ ,(4)»k^k-w\ ,»,, , 

^int.O — | 4' 5 f \ int >° int.O ^int.O ^int.O J f l K ) w J) 

C% = - [\, I] (42,0 + 4nto) (V (k',w)0* (k',q) - 0*(k' - q,w)^(k' - w,q)), 

^2,0 =~{i\} Kt )k o'^ + V ( k ' + U - * U ) - ^ + U, U)) , 

<o = - { ^ } (<o^ (k - w, u) - k " V (k - q, u)) , 

*ffio = - { \\ } Sg (0* (1, u)0(l, w) - ^ (1 - w, u)^(l - u, w)) , 
G&, = {^} Kio^-Xk' Q, ») BS2# «)) , 
*&> = {\,\} «r k "^(k - w, u) - fl« #k - q, u)) , 
4go = {\>l} B So (^ + w > "W 1 - u + w, w) - 0(1, u)0(l + w, w)) , 
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K 



(5) 
int.O 



1 1 

4' 5 



Cg ^(k - q, u) - Cggr*" V(k + w, u)) 



-(5) 
J int,0 



{i' J} ( C So^(k',u) - c£ k o'^ k '->*(k' - q,u)) 



M (5) - 

iW int,0 — 



1 1 

4'5 



Cg (f (1, u)0(l, w) - <P* (1 - w, u)0(l - u, w)) , 



(5) 



int,0 



J. ^} (^Sr k -^(k - q, u) - Cffiotfk + w, u)) , 



©So = {^} (Cgo^ +«,«)- CW^W + u - q, u)) 



p(5) . 
Mnt.O _ 1 4 



J, 1 1 Cg + w, u)0(l - u + w, w) - 0(1, u)0(l + w, w)) , 



^int.O 1 4 ' 5 f V int '° int .° / ' °int,0 



1 1 



1 1 K E {i) T (5) 

4' 5 | int.O' -MntjO 



4' 5 J V mt, ° mt >° / 



(5) 



1 1 



^(4) . r ,(4)*k^k' 



(•>) 



'int,0 _ "I 4' 5 f V int >° int,0 J ' K, 'mt,0 _ 1 4 

1/4 in the curly brackets refering to int terms and 1/5 to terms 



III (/.W +/ (4)*k^k'\ 
4' 5 I \ int '° in t.O y ' 



3. Sixth order 



H(6) = l4 [5 ' ^ int]]]]] + 7^0 [5 ' [5 ' [5 ' [S ' [S > [S ' H ° 



— ^int + ^0 . 

(C4) 
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rr(6) 



o = { g. £ } | H AS,o7w7q«ka k '_ w ak'+q + XI S int,o7w7qakak-w7q + Qnt,o7 wak '- wak '-q nk 7q 

I qkwk' qkw qkwk' 

i ^int,o7u7wa k _ w ak- q a k /_ q a k '+u + Ant,o7u7wa k _ w _ u a k _ q a k /_ q a k > 

qkwk'u qkwk'u 

Y F intV ua l-« ai - wa k-w a k-q Cl k'-q a k' + ^ G inlo7wa k _ w a k _ q a k ,_ q a k ,_ u 7* 
qkwk'ul qkwk'u 

/ i„t,o7wa k _ w a k _ q _ u a k ,_ q a k - 7u + ^ •AntVi+w a i-ua k _ w a k _ q a k ,_ q a k ,7 u 

qkwk'u qkwk'ul 

Y ^int,o7u7wa k ,a k -_ q a k _ q _ u a k+w + ^ o7u7w a k ,_ u a k /_ q a k _ q a k+w 

qkwk'u qkwk'u 

X! M int ) ,o7uaJ'_ u ai_ q a k ,a k /_ q a k _ q a k , +w + ^ iV i ( n 5 t ) 7 w a k ,a k /_ q a k _ q a k+w _ u 7* 

qkwk' ul qkwk' u 

X! O i ( nt,o7wa k '+uak'-qa k _ q a k+w 7* + ^ P i ( I f t ) a 1 * +w ai_ u a k ,a k /_ q a k _ q a k+w 7* + ^ ^t,o7q"-k«k'7 q 

qkwk'u qkwk'ul qkk' 

' S int,07qak a k-q7 q 7 q + 2^ <S mt,07qa k ak-q7 q 7 q + 2^ T int,o7wak-w tt k-qffl k '-qOk'-w7w 



$1 C/ intV wa k'+w a k'-qa k _ q ak+w7w + Y M/ int ) ,o7wa k a k _ q a k ,_ q a k '7^ + 6*6{} + c, 

qkwk' qkwk' 



(C5) 



The last term denoting all (6, 2) expressions which are irrelevant after averaging over phonon vacuum. The coefficients 
are the same as in the 5th order terms, except for the factors 1/5, 1/6. 
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